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Abstract We used recently produced Solar System ephemerides, which incorporate two 
years of ranging observations to the MESSENGER spacecraft, to extract the secular orbital 
elements for Mercury and associated uncertainties. As Mercury is in a stable 3:2 spin-orbit 
resonance these values constitute an important reference for the planet’s measured rotational 
parameters, which in turn strongly bear on physical interpretation of Mercury’s interior 
structure. In particular, we derive a mean orbital period of (87.96934962 ±0.00000037) days 
and (assuming a perfect resonance) a spin rate of (6.138506839 ±0.000000028)°/day. The 
difference between this rotation rate and the currently adopted rotation rate (Archinal et al, 
2011) corresponds to a longitudinal displacement of approx. 67 m per year at the equator. 
Moreover, we present a basic approach for the calculation of the orientation of the instan¬ 
taneous Laplace and Cassini planes of Mercury. The analysis allows us to assess the uncer¬ 
tainties in physical parameters of the planet, when derived from observations of Mercury’s 
rotation. 

Keywords Mercury ■ spin-orbit coupling ■ Laplace plane • MESSENGER • ephemeris 
PACS 96.30.Dz • 96.12.De ■ 95.10.Km 


1 Introduction 

Mercury’s orbit is not inertially stable but exposed to various perturbations which over long 
time scales lead to a chaotic motion (Laskar, 1989). The short-term (about few thousand 
years) evolution of the orbit can be approximated by a secular contribution to the orbital ele¬ 
ments. Most prominent is the precession of the pericenter of Mercury’s orbit, which was also 
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Fig. 1 The unit vectors ex.y.z denote the International Celestial Reference Frame (ICRF). Mercury's orbital 
plane is illustrated by a dashed ellipse and its orientation by the vector e 0 . The ecliptic and the ex _e Y plane 
of the ICRF are given by dotted and dash-dotted ellipses, respectively. The Laplace plane normal is indicated 
by e„ and Mercury's spin axis by e s . The figure is based on numbers given in Tab. 1 at the J2000.0 epoch. 


an important test of Einstein’s theory of general relativity (Einstein, 1915). Due to the Sun’s 
torque on the asymmetric mass distribution of Mercury, the rotation of Mercury is strongly 
coupled to its evolving orbit. Radar observations (Pettengill and Dyce, 1965) revealed that 
Mercury’s rotation period is about 59 days and in a stable 3:2 resonance with its orbital 
period (Peale and Gold, 1965; Colombo, 1965). More recently Margot et al (2007) have 
used an Earth-based radar-speckle correlation technique to precisely measure the physical 
libration amplitude and the obliquity of Mercury. By interpretation of these measurements 
in terms of physical parameters of the planet - following the idea of the Peale experiment 
(Peale, 1976, 1981) - the authors concluded that Mercury’s core is at least partially molten 
(Margot et al, 2007, 2012). 

With the MESSENGER space probe (MErcury Surface, Space ENvironment, GEo- 
chemistry, and Ranging) having entered orbit around Mercury in March 2011, the obser¬ 
vational data of Mercury have greatly improved. Further, new Solar System ephemerides 
which incorporate two years of ranging and Doppler tracking observations to MESSENGER 
were produced. For the interpretation of the observations of Mercury’s rotation performed 
by instruments on MESSENGER, precise knowledge of the resonant rotation parameters of 
Mercury is mandatory. In fact, the resonant spin rate, currently adopted in the rotation model 
of Mercury, dates back to the first IAU report (Davies et al, 1980). 

In this work we provide updated reference values for Mercury’s rotation assuming the 
perfectly resonant rotation model based on the most recent planetary ephemerides. These 
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values serve as a basis for the interpretation of the rotational parameters of Mercury, which 
are proposed to be measured with high precision (Stark et al, 2015). 


2 Secular orbital elements of Mercury 

Recently, new Solar System ephemerides DE432 from the Jet Propulsion Laboratory (W. 
M. Folkner, personal communication, 2014) and INPOP13c from the Institut de Mecanique 
Celeste et de Calcul des Ephemerides (Fienga et al, 2014) were produced. Besides other 
improvements these ephemerides incorporate updates to the orbit of Mercury. Both ephe¬ 
merides although different in their production process and covered time span led to identical 
results in our calculations. We concentrate here on the DE432 ephemeris and give the orbital 
elements derived from the INPOP 13c ephemeris in appendix 4. 

The DE432 ephemeris covers a time span of approximately 1,000 years (1 January 1550 
to 1 January 2550). In this time span we derived the osculating Keplerian orbital elements 
of Mercury from state vectors given with respect to the Sun-centered International Celestial 
Reference Frame (ICRF). We used a time step of 7 days and set the gravitational parameter 
of the Sun to GMq = 132,712,440,041.9394km 3 /s 2 (Folkner et al, 2014). For the calcu¬ 
lation of the osculating orbital elements standard techniques were used (Bate et al, 1970). 
In order to obtain the secular parts of the elements we decomposed the osculating orbital 
elements in a quadratic polynomial and a sum of periodic terms 

x{t ) =xo+xit/cy + x 2 {t/cy) 2 + Y J AiCos(v i t + (l>i), (1) 

i 

where x stands for a Keplerian orbital element a, e, /, Q , ft), M, being semi-major axis, eccen¬ 
tricity, inclination, longitude of ascending node, argument of pericenter, and mean anomaly, 
respectively. The time t is measured in Julian centuries (cy) from the J2000.0 epoch. Higher 
order terms in the polynomial were discarded because their estimated uncertainty exceeded 
the actual value by many orders of magnitude. The periodic terms are characterized by their 
amplitude Aj, frequency v,, and phase </),. We list these values for ten highest amplitudes of 
each orbital element in appendix 4. 

The decomposition of the osculating orbital elements time series into the form of Eq. 
1 was performed with the help of the frequency mapping tool FAMOUS 1 . This is done 
because a simple least-squares fit may lead to biased results given the fact that the variations 
of the orbital elements are in first order periodic and not random. At least 50 frequencies 
were identified and subtracted from the variation of the orbital elements. The variance of the 
periodic variations <7 2 was used to estimate an uncertainty for the coefficients. Thus, orbital 
elements with relatively high periodic variations receive higher error bars. The uncertainties 
of the secular coefficients x\ and x 2 were derived by considering the maximal slope and 
curvature of the polynomial within the interval [—(Tv, <tJ and a time span of 1000 years. The 
resulting values are given in Table 1. 

In order to demonstrate the convergence of the method we increased the number of 
frequencies to 100 and found only changes below 2% of the uncertainties of the polyno¬ 
mial coefficients. For further verification of our approach we calculated orbital elements 
with respect to the ecliptic at J2000.0 (see appendix 4) and compared our results with those 
published by Standish and Williams (2013). Beside the secular parts of the inclination and 
longitude of ascending node our values and their uncertainties are consistent with the pub¬ 
lished values. The discrepancy we found in I\ and Q\ may result from the fact that we 

1 F. Mignard, OCA/CNRS, ftp://ftp.obs-nice.fr/pub/mignard/Famous 
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x 0 

Xi 


a/(10 6 km) 

57.90909 

0.002 x 10~ 6 

-0.002 x 10 


±0.00011 

±22.34 x 10~ 6 

±4.45 x 10 

e 

0.2056317 

20.4 x 10~ 6 

-20 x 10 


±0.0000071 

±1.4 x 10~ 6 

±290 x 10 

i/° 

28.552197 

0.0048464 

-9.8 x 10 


±0.000036 

±0.0000073 

±1.5 x 10 

n/° 

10.987971 

-0.032808 

-12.3 x 10 


±0.000099 

±0.000020 

±4.0 x 10 

(o/° 

67.5642 

0.18861 

-3 x 10 


±0.0020 

±0.00040 

±80 x 10 

M/° 

174.7948 

149472.51579 

8x 10 


±0.0032 

±0.00063 

±126 x 10 


Table 1 Secular Keplerian orbital elements of Mercury as derived from the DE432 ephemeris at epoch 
J2000.0, given with respect to ICRF (see Fig. 1). 


consider the quadratic term, which is significant for these elements. Comparison with other 
literature values (Margot, 2009; Noyelles and D’Hoedt, 2012; Noyelles and Lhotka, 2013) 
shows excellent agreement with our values for these orbital elements. 

Additionally, we calculated the precession of the pericenter of Mercury. Note that the 
secular rates are strongly dependent on the selected reference frame. We used the mean 
orbital plane of Mercury at J2000.0 (see Sec. 2.2) as reference frame and found a precession 
of 575.3 ± 1.5arc sec/cy (see appendix 4). Again this is in a very good agreement to the 
literature value of (5600.73 — 5025) arc sec/cy = 575.73 ± 0.41 arc sec/cy (computed from 
Weinberg, 1972, p.199). 

Another method to obtain the secular orbital elements involve the usage of a secular po¬ 
tential and integration of the averaged differential equation of Mercury’s motion. However, 
such a method neglects the mutual interaction of the perturbing planets and is not appropri¬ 
ate for precise interpretation of spacecraft data (Yseboodt and Margot, 2006). More details 
on averaging methods can be found in e.g., Sanders et al (2007). 


2.1 Mean orbital period 

The mean period of the orbit is defined as ^orbit = Itz/uq, where no is the mean motion of 
Mercury. We can derive no from the first order term of the mean anomaly M = Mo +M\t = 
no(to + t). The time to is the elapsed time at J2000.0 since the last pericenter passage. Using 
the values in Tab. 1 we derive 


n 0 = M x = (4.092334450 ± 0.000000017)°/day (2) 

to = Mq /M, = (42.71274 ± 0.00077) day (3) 

7’orbit = 360°/Mi = (87.96934962 ± 0.00000037) day. (4) 


In order to check the derived value of no we used Kepler’s third law no = JGMq/oq and 


found 4.092343 ± 0.000083 °/day. This value is consistent with Eq. 2 but has an error larger 
by two orders of magnitude. 
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2.2 Mean orbital plane 

From the secular parts of the orbital elements Q and I we can calculate the mean normal 
vector of the orbital plane by 

e 0 = sinTl sin/ex —cosX2 sin/e Y ±cos/ez (5) 

= cos a orbit cos S orbit ex + sin a orbit cos <5 orbit e Y + sin S orbit e z , (6) 

where ex.y.z denote the orientation of the ICRF (see also Fig. 1). Comparing Eq. 5 and 
6 we find the right ascension and declination of the orbit pole to be a olbb = Q — tz/2 
and <5 orbit = n/2-L At J2000.0 the values are a° rbit = (280.987971 ± 0.000099)° and 
<5° lbb = (61.447803 ± 0.000036)°, respectively. From the secular components of Q. and 
I we can directly derive the first order precession rates of the orbit pole a \ )lbb — Q\ = 
(—0.032808±0.000020)°/cy and S° rbit = ~h = (-0.0048464±0.0000073)°/cy. It should 
be noted that the precession of the orbit normal is treated here as a secular variation in 
inclination and longitude of ascending node, which is justified by the long period of the 
precession. By that reason the given description of the mean orbital plane is strictly valid 
only for the time span of the ephemeris, i.e. about ±500 years around the J2000.0 epoch. 
The error bars on the orbit pole orientation and precession rates were obtained through 
propagation of the uncertainties in the orbital elements. Note that the derived values are in 
agreement with the findings of Margot (2009) with a orblt = 280.9880° — 0.0328°t/cy and 
gorbit _ 61 .4478° — 0.0049°t/cy, where DE408 ephemeris and a period of 200 year was 
used. 


2.3 Laplace plane 

The other planets of the Solar System exert a torque on the orbital plane of Mercury, which 
leads to a quasi-periodic precession of the orbit normal. Further, the plane to which the incli¬ 
nation of Mercury remains constant, i.e., the Laplace plane, also undergoes slow variations 
(Noyelles and D’Hoedt, 2012). Several attempts have been made to calculate the orienta¬ 
tion of the Laplace plane normal (Yseboodt and Margot, 2006; Peale, 2006; D’Hoedt et al, 
2009), each of them leading to different results in the Laplace pole position and the preces¬ 
sion period (see Fig. 2). 

The concept of the ’’instantaneous” Laplace plane was proposed for Mercury by Yse¬ 
boodt and Margot (2006) to derive an approximate Laplace plane valid for several thousand 
years. Note that without additional assumptions the instantaneous precession vector w is 
only constrained to a line. In order to overcome the ambiguity in the instantaneous Laplace 
plane either a fit to the ephemeris (Yseboodt and Margot, 2006) or some additional assump¬ 
tions (Peale, 2006; D'Hoedt et al, 2009) are used. 

Here we introduce a concept of an instantaneous Laplace plane, which removes the am¬ 
biguity in its instantaneous orientation and precession frequency. As the only assumption 
we require the instantaneous Laplace plane to be invariable, i.e. w = 0. Note that a similar 
concept was suggested by Yseboodt (2011). When an instantaneous Laplace plane is consid¬ 
ered it is also important to clarify in which time period it should be instantaneous. The free 
precession period of Mercury is in the order of 1000 years (Peale, 2005). This means that 
the rotation axis will not be affected by the short period (on the order of decades) changes 
in the orientation of the orbit normal, but will follow the changes at long periods due to 
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adiabatic invariance (Peale, 2005). In order to obtain an instantaneous Laplace plane, which 
is relevant for Mercury’s spin, we can neglect all periodic variations and consider only the 
secular terms. 

The general equation for the precession around an axis w is 

w x e„ = e„. (7) 


The precession vector w is given by w = —fi e w , where e w is the orientation of the Laplace 
plane and /I the precession rate. First we multiply both sides of Eq. 7 with e 0 and obtain 

w = (e 0 x e 0 ) - At cost e 0 , (8) 


where we used e 0 ■ w = — At cos I and t is the inclination of Mercury’s orbit with respect to the 
Laplace plane. In order to constrain the instantaneous orientation of the Laplace plane we 
have to find an instantaneous value for fi cost. By differentiating Eq. 7 and the requirement 
w = 0 we obtain 

e 0 + F 2 e o = —jttcosiw, (9) 

where we make use of w x e 0 = w x (w x e 0 ) = —At cos I w — H 2 e 0 - The differential equation 
Eq. 9 describes a regular (uniform) rotation of e 0 around e w with the frequency At. Following 
the concept of the instantaneous Laplace plane Eq. 9 is only fulfilled with a unique w for a 
specific time t. By multiplying Eq. 9 with e 0 and using Eq. 7 we can find 


At cosi = 


e„ ■ (e„ x e„) 


( 10 ) 


By using Eq. 5 for e 0 and e 0 we can obtain the values of n cos I as it would be in a regular 
form and by that an expression for the instantaneous Laplace plane. The combination of Eq. 
8 and 10 gives the instantaneous Laplace plane orientation 


w = (e 0 x e„) 


e„ ■ (Go x Ho) 


( 11 ) 


Following the formalism of Peale (2006) Eq. 11 can be expressed as 


w= (/cos/2 + (w z — Ll) tan/sin/2) ex + 

+ (/sin/2 — (w z — /2) tan/cos/2) ey + w ; ez (12) 


and w z given by 


w z — /2 -t- 


(IQ — //2)sin/ + /2/ 2 cos/ 

-r;-■-r-cos / 

I 2 + (Z2 sin/) 2 


(13) 


The instantaneous Laplace pole given by Eq. 11 is practically equivalent to the fit of the 
ephemeris to a cone, performed by Yseboodt and Margot (2006). 

Using Eq. 11 we calculate the coordinates of the Laplace pole at J2000.0 to a^ p = 
(273.8± 1.0)° and 5 PP = (69.50±0.77)° in the ICRF. It should be noted that the covariance 
Cov(a pp , <5 PP ) = —(0.77°) 2 is very high, indicating a high correlation (—99.8%) of the 
right ascension and declination values. The instantaneous precession rate is At = (0.00192 ± 
0.00018)/cy, 7 lp = (327300± 32000)years and i = (8.58±0.84)°. 

However, if Mercury is in a Cassini state (see Sec. 3.4) the determination of the po¬ 
lar moment of inertia from the obliquity is not largely affected by the uncertainties of the 
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Fig. 2 Orientation of the orbit (OP) and the Laplace plane (LP, black disk with error bars) normal at J2000.0 
epoch with respect to the ICRF. The precession of the orbit pole around the instantaneous Laplace pole is 
indicated by a dashed arc. The Laplace plane orientation of Yseboodt and Margot (2006) (grey disk) overlaps 
with the values derived in this work. Note that the longitude of the Laplace pole is given incorrect in (D'Hoedt 
et al, 2009) and was corrected in (Noyelles and D’Hoedt, 2012). The figure shows the corrected position (blue 
disk). 


Laplace pole (Yseboodt and Margot, 2006; Peale, 2006). Using our formalism we can derive 
the relevant quantities and their errors 

/isint = (2.8645±0.0016) x 10~ 6 /years (14) 

jUcosi= (18.98 =L 1.83) x 10~ 6 /years. (15) 

The correlation between jiisini and jUcost is very low Corr(jU sint,/t cost) = —10~ 3 . 


3 Mercury’s rotation 

3.1 Rotation model 

The rotation model of a celestial body consists of a set of values defining its orientation as a 
function of time with respect to a reference frame. Here we recall briefly the IAU convention 
of a rotation model for Mercury (e.g. Archinal et al, 201 1). 

The orientation of Mercury’s spin axis is described by the right ascension a and decli¬ 
nation 5 coordinates of the intercept of the spin axis vector e s with the celestial sphere. The 
spin axis vector e s with respect to the ICRF is given by 


e s = cos a cos 5 ex + sin a cos 5 ey + sin 5 ez, 


(16) 
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where ex.Y.z denotes the ICRF. The orientation at J2000.0 epoch is denoted by do and So 
and the first order precession rates are (X\ and d|. The rotational axis is consequently given 
by 

a(t) = ao + a\t (17) 

5(r) = 5 0 + 5 I t. (18) 

The rotation of Mercury around its axis is described by the longitude of the prime meridian 
<po. the rotation rate <pi, and the physical longitudinal libration <jpm, 

<P( f ) = < Po + <P 1 t + <Pi ib (0- (19) 

The rotation model is of great importance, as it is used to derive body-fixed coordinates 
of observations performed by spacecraft. The matrix R which transforms coordinates from 
ICRF to body-fixed is composed from three rotations 

R = Rj(p)Rj(jr/2-5)Rj(a + jr/2), (20) 

where R x , z denote counter-clockwise rotations (right hand rule) around the x- and z-axis, re¬ 
spectively. Note that the spin axis orientation can be computed by e s = R z (a±7r/2)R x (7r/2 — 
<5)ez- 


3.2 Resonant rotation 

Using the secular orbital elements in Tab. 1 we derive the resonant spin rate of Mercury 

13 12 ) 

<pf for the case that the spin is in perfect 3:2 resonance to the motion of the planet on 
its orbit. Using the mean motion value no = M\ derived from the mean anomaly and the 
precession of the argument of pericenter ft>i we compute the mean resonant spin rate to 

(p{ 3/2) = ^no + coi = (6.138506839 ± 0.000000028) 0 /day. (21) 

The current value of Mercury’s rotation found in the literature is 6.1385025°/day (Archinal 
et al, 201 1). The difference between these rates corresponds to a longitudinal displacement 
of 5.7arc sec per year (approx. 67m per year at the equator of Mercury), which should be 
noticeable during e.g. the MESSENGER mapping mission (where typical image resolution 
vary from few kilometers to few meters). 

(3 / 2 ) 

We want to stress that the resonant spin rate <p) ' in the rotation model is composed 
of the planet’s rotation around its spin axis and the precession of Mercury’s orbit. Thereby, 
we have to consider the precession of the argument of pericenter 0 >i = (5.164 ±0.011) x 
10 -6o /day, and not the precession of the longitude of pericenter GJi = ± ®i since the 

precession of the ascending node Q\ is already incorporated in the precession of the rota- 

13 / 2 ) 

tional axis. Note that the spin rate (p) ' Ms defined with respect to a precessing frame and is 
not strictly ’’sidereal” since the rotation axis changes slowly its orientation. Further, Mercury 
has a small but non-zero obliquity of i c = 2.04arc min (Margot et al, 2012). However, the 
correction arising from the obliquity is in the order of i c Cl^/I\ and can be neglected when 
comparing with the error of ft>i (see appendix 4). 

If one of the sub-solar points at perihelion is used for the definition of the prime meridian 
<po, the orientation of Mercury’s long axis at J2000.0 with resonant rotation would be 

(p£ /2) = ^M 0 + coo = ( 329 . 7564 ± 0 . 0051 )°. 


(22) 
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We find excellent agreement of this value with the findings of Margot (2009), who stated a 
value of 329.75°. Note that the actual prime meridian of Mercury tpo is defined with respect 
to the crater Hun Kal located at 20° W (Archinal et al, 2011). 


3.3 Physical libration in longitude 


The annual libration of Mercury is closely tied to the revolution of Mercury around the Sun. 
Of particular importance is the mean anomaly M of Mercury since it defines the period and 
the phase of the libration. The libration is modeled as follows (Goldreich and Peale, 1966) 


9lib( f ) =£g88/ J fSin(£no(f + to)) • 

k 

The amplitudes follow a recursive relationship 

G2oi(£+ 1, eo) 

*88/(* +1 )-*88/* G20l(Mo) - 

where G 2 oi{k,eo) are given by Kaula’s eccentricity functions (Kaula, 2000) 


G20\{k,e) = 


G20\-k{e) — G2oi+*r(e) 
k 2 


Using eo = 0.2056317 ±0.0000071 from Tab. 1 we calculate the first five terms to 


(23) 


(24) 


(25) 


G 2 oi(l,eo)= 0.569650 ±0.000027 (26) 

G 2 oi (2,£ 0 ) = (-60.0733 ±0.0042) x 10~ 3 (27) 

G 201 (3,£ 0 ) = (-5920.32 ±0.77) x 10~ 6 (28) 

G 201 (4,£ 0 ) = (-1200.10 ±0.20) x 10“ 6 (29) 

G 201 (5,£q) = (-267.691 ±0.053) x 10~ 6 . (30) 


The main period of the annual libration is the mean orbital period 7),,.^ (Eq. 4) with 
the phase given by Mq = ;?o/o- In addition, long-period variations of the orbital elements can 
lead to forced librational motion of Mercury with periods other than the orbital period (Peale 
et al, 2007; Yseboodt et al, 2010), but these are not considered in this work. 

The measurement of the libration amplitude provides important constraints on the in¬ 
terior structure of Mercury. The amplitude of the annual libration g88 is related by (Peale, 
1981) 

g&& = k ) —^ —G201 (l±o) ( 31 ) 

z c m 

to the ratio of moments of inertia (B — A)/C m , where A < B < C are the principal axes of 
inertia of the planet and C m is the polar moment of inertia of the mantle and crust. Assuming 
the libration amplitude could be measured with a negligible error, the uncertainty in (B — 
A)/C m would be only at 6 x 10~ 7 , due to the uncertainty in the eccentricity of Mercury’s 
orbit. Here we used a libration amplitude of gs8 = 38.5 arc sec (Margot et al, 2012). 
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3.4 Cassini state 


Mercury is assumed to occupy a Cassini state 1 (Peale, 1969), implying that the spin vector 
of Mercury e s lies in the plane defined by the Laplace plane normal e w and the orbit normal 
e 0 with the latter being enclosed by the others. The spin axis is consequently in a 1:1 reso¬ 
nance to the precession of the orbit normal, i.e. CC\ ~ a” 1 ” 1 and 5i « S”* 3 ' 1 . Note that the 
spin axis precesses with slightly higher rates as described in appendix 4. The Cassini plane 
e c , which contains all Cassini states can be expressed as a linear combination of the orbit 
and Laplace plane normal 

e c = re 0 + .sw. (32) 

We can constrain r and 5 by je c | = 1 and e c ■ e 0 = cos i c , where i c is the obliquity. By using 
Eq. 8 for w and je 0 x e 0 | = |e 0 | this results in 


with 


, . . e 0 x e„ 
e c = cos i c e 0 + sin i c ———— 

|e 0 1 


(33) 


e 0 x e 0 = (/cos/2 — Q sin/cos / sin /2) ex 

+ (/sin/2 + /2 sin/cos/cos/2) ey + /2(sin/) 2 ez (34) 

e 0 | = jUsint = \Ji 2 + (Q sin/) 2 . (35) 

From Eq. 33 it can be verified that the plane defining the Cassini state is independent 
from the exact form of the precession of the orbit around the Laplace plane. Especially, it is 
not dependent on w z as recognized by Peale (2006) and Yseboodt and Margot (2006). The 
Cassini plane is sufficiently defined by the orientation of the orbit normal and its temporal 
change. In fact, the Cassini plane normal is the vector e 0 given by 

e 0 = (/2 cos/2 sin/ + /sin/2 cos/) ex 

+ (/2sin/2 sin/ — /cos Q cos/) ey (36) 

— /sin/ez ■ 


With the uncertainty of the orbital elements we can estimate the ’’thickness” of the 
Cassini plane, which results from uncertainties in the knowledge of the secular variation 
of Mercury’s ephemeris. At i c = 2.04arc min (Margot et al, 2012) we find a lc thickness of 
0.18 arc sec. This allows to interpret possible offsets of Mercury’s spin orientation from the 
exact Cassini state (Margot et al, 2012; Peale et al, 2014). 

Using the obliquity i c the polar moment of inertia C/ mR 2 (scaled with the mass m and 
radius R of Mercury) can be calculated by (Peale, 1981) 

C _ n 0 sini c ((7 2 (l -e 2 )~ 3//2 cosi c +C 22 G 2 Qi(g)(l +cosi c )) 
mR 2 jii sin t cos i c — /.I cos t sin i c 

where Ji = 5.03216 x 10~ 5 and C 22 = 0.80389 x 10~ 5 (Mazarico et al, 2014) are the second 
degree harmonic coefficients of Mercury’s gravity field. Assuming perfect knowledge in the 
obliquity and the gravitational coefficients the error on C/mR 2 is only at 6.1 x 10 5 due 
to the uncertainty of the orbital elements. Note that our analysis does not include model 
uncertainties of Eq. 37, e.g., simplifying assumptions which were made in the derivation of 
the equation. We only infer the uncertainty of C/mR 2 due to orbital elements if Eq. 37 holds 
exactly and all other quantities are perfectly known. A more sophisticated analysis including 
higher order gravity field and tides can be found in Noyelles and Lhotka (2013). 
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4 Discussion and conclusion 

In this work we extract orbital elements for Mercury from ephemeris data and predict a mean 
rotational model for Mercury in the view of a perfect resonance to its orbit. In this case the 
rotation is, besides the obliquity and the libration amplitude, completely determined by the 
mean orbital elements and their rates. On the basis of the uncertainties in the mean orbital 
elements, errors of the theoretical perfectly resonant rotation model can be estimated. Note 
that ephemeris uncertainties are estimated from the periodic variation of the orbital elements 
and do not reflect any accuracy or ’’error” of the ephemeris. They can be rather understood as 
model uncertainties, since the secular part of Mercury’s orbital elements does not capture the 
full variation of the orbit. In this work we introduced a consistent approach which allows us 
to estimate the uncertainties of the rotational parameters resulting from the simplified secular 
orbital elements. The findings are of great importance for interpretation of the current and 
future observations of Mercury’s rotation by MESSENGER and BepiColombo spacecraft. 


Appendix 1 

The Keplerian orbital elements derived from the INPOP 13c ephemeris (Fienga et al, 2014) 
are given in Tab. 2. We find very little difference of the values when comparing to the DE432 
ephemeris (see Tab. 1). The deviation for the trend a\ of the semi-major axis is about 1.7 
meter per century. 



xo 

XI 

*2 

a/(10 6 km) 

57.90909 

—44 x 10~ 12 

420 x 10- 12 

e 

0.2056317 

20.4 x 10~ 6 

-0.027 x 10~ 6 

i/° 

28.552197 

0.0048473 

-9.8 x 10~ 6 

n/° 

10.987969 

-0.032808 

-11.9 x 10~ 6 

(D/° 

67.5642 

0.18862 

—4 x 10~ 6 

m/° 

174.7948 

149472.51578 

7x 10- 6 


Table 2 The same as Tab. 1 but derived from the INPOP13c ephemeris and a time span of 2000 years 
(09.06.973 AD - 23.06.2973). 


In Tab. 3 we give values for reference frame dependent orbital elements with respect 
to the ecliptic (ECLIP, inclination 23.439291°), Mercury orbital plane (OP), and Mercury 
Laplace plane (LP) at J2000.0. The rotation matrices for the transformation to the these 
reference frames from the ICRF are given by 


1 0 0 \ 

R eclip = 0 0.91748206 0.39777716 , 

\0 -0.39777716 0.91748206 / 

/ 0.98166722 0.19060290 0\ 

R op = -0.16742216 0.86227887 0.47795918 , 

\ 0.09110040 -0.46919686 0.87838205 ) 

/ 0.88845611 0.43672271 0.14113473 \ 

R lp = -0.45838720 0.82896828 0.32045711 . 

\ 0.02295468 -0.34940643 0.93669004 / 


(38) 


(39) 


( 40 ) 
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The precession of the pericenter of Mercury is given by = 72° p + 0 >° p = 575.3 arc sec/cy. 

The inclination of the orbital plane with respect to the Laplace plane 7 q p = t = 8.58° remains 
constant 7[ p = Pf ~ 0. The precession of the orbit around the Laplace plane is 
0.109981°/cy. ' 


m=n = 



XQ 

X\ 

*2 

7/° 

7.004975 

Ecliptic 

0.0059524 

0.7 x 10~ 6 

72/° 

48.330908 

-0.125416 

-89.2 x 10~ 6 

co/° 

29.1252 

0.28428 

80 x10~ 6 

m/° 

77.4561 

0.15886 

-13 x 10~ 6 

l/° 

Mercury orbital plane 

0.0 -0.016413 -3.9 x nr 6 

n/° 

68.735669 

-0.054375 

337.0 x 10~ 6 

co/° 

320.3895 

0.21417 

-350x10~ 6 

®/° 

29.1252 

0.15980 

-13 x 10~ 6 

i/° 

Mercury Laplace plane 

8.582338 1 x 10~ 18 5 x 10“ 21 

n/° 

0.0 

-0.109981 

-25.9 x 10~ 6 

<o/° 

50.3895 

0.26855 

13x10~ 6 

m/° 

50.3895 

0.15857 

-13 x 10~ 6 


Table 3 Orbital elements of Mercury as derived from the DE432 ephemeris at epoch J2000.0 with respect to 
the following reference frames: Ecliptic and Earth equinox of J2000; Mercury orbital plane of J2000.0 and 
ascending node with respect to the ecliptic; Mercury Laplace plane and ascending node with respect to the 
Mercury orbital plane of J2000.0. 


Appendix 2 

In Tab. 4 we list the first ten periodic terms, which were identified in the osculating orbital 
elements time series. Some of the periods can be assigned to planetary perturbations, e.g., 
Venus: (Ay) 0.62 years; (2Ay) 0.31 years; (2A m — 5Ay) 5.66 years; (A m — 3Ay) 1.38 years; 
(Am —2Ay) 1.11 years; (2Am — 4Ay) 0.55 years; Earth: (Am — 4A e ) 6.58 years; Jupiter: 
(Ay) 11.86 years; (2Ay) 5.93 years; (3Ay) 3.95 years; Saturn: (2A^) 14.73 years, where 
A = M + GJ = M + Q + co denotes the mean longitude of the planet, respectively. 


Appendix 3 

The obliquity of the spin axis i c introduces small changes in the precession and resonant 
rotation rates. To stay within the Cassini plane the spin axis has to precess slightly faster 
than the orbital plane normal. In order to compute the corrections we expand equation Eq. 
33 to first order in the obliquity i c . The declination 5 and right ascension a of the spin axis 
are then given by 

, n 72 sin7 . _ 

5(f) — — — 7H— i c — 5o + 5if 

y 7 2 + (72 sin7) 2 


( 41 ) 
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i 

a [km] 

A,- 

T, 


e [10 6 
A,- 

] 

T, 

<t>i 

1 

109.56 

1.11 

145.25 

7.23 

5.93 

272.97 

2 

56.37 

0.20 

140.65 

4.63 

1.11 

326.10 

3 

54.55 

5.66 

356.60 

3.53 

5.66 

181.05 

4 

35.31 

0.29 

213.41 

2.53 

1.38 

141.21 

5 

31.01 

0.12 

76.038 

1.52 

11.86 

123.83 

6 

29.41 

1.38 

325.51 

0.89 

0.25 

261.39 

7 

21.92 

0.13 

211.50 

0.88 

0.55 

109.00 

8 

21.19 

0.25 

81.45 

0.86 

14.73 

304.98 

9 

20.27 

0.55 

285.27 

0.81 

0.46 

128.90 

10 

19.46 

0.40 

70.013 

0.78 

0.29 

36.27 


/ [10 3 arc sec] 


n [io 

3 arc sec] 


i 

A; 

T t 

<t>i 

A; 

Ti 

<t>i 

1 

167.3 

5.93 

15.01 

399.8 

5.93 

292.75 

2 

52.5 

5.66 

71.86 

165.8 

5.66 

135.33 

3 

31.9 

1.38 

250.97 

145.2 

11.86 

145.14 

4 

22.1 

11.86 

267.19 

116.6 

1.11 

249.76 

5 

20.1 

14.73 

58.16 

107.3 

1.38 

189.18 

6 

17.5 

6.58 

343.74 

57.5 

0.62 

356.75 

7 

17.2 

3.95 

51.01 

52.3 

0.40 

155.89 

8 

11.8 

0.31 

326.19 

48.3 

14.73 

336.52 

9 

9.9 

0.24 

112.99 

44.0 

6.56 

59.20 

10 

9.6 

0.12 

99.47 

42.6 

0.24 

45.75 


co [arc sec] 


M [arc sec] 


i 

A/ 

T, 

1 Pi 

A/ 

T, 

<l>i 

1 

7.36 

5.93 

180.64 

10.71 

5.66 

87.97 

2 

4.57 

1.11 

55.93 

8.04 

1.11 

235.40 

3 

3.49 

5.66 

272.57 

7.70 

5.93 

3.50 

4 

2.55 

1.38 

50.09 

1.92 

1.38 

230.24 

5 

1.62 

11.86 

17.21 

1.90 

11.86 

186.46 

6 

0.89 

14.73 

212.80 

1.31 

6.57 

334.16 

7 

0.84 

0.25 

351.58 

1.22 

0.55 

17.86 

8 

0.84 

0.55 

200.04 

1.10 

0.46 

38.25 

9 

0.78 

0.46 

219.04 

1.10 

0.29 

305.42 

10 

0.78 

0.62 

194.46 

1.09 

0.25 

171.25 


Table 4 Ten leading terms of the decomposition of the time series of the osculating orbital elements in 
£, AjCOs(Vjf + 0/) with v,- =2iz/Ti. The unit of the amplitude A,- is given in the brackets beside each orbital 
elements symbol, respectively. The periods 7} are given in years and the phases 0, in degrees. The values are 
given for orbital elements in the ICRF and derived from the DE432 ephemeris. 


and 


a(t) = £2-- + 


7/ sin/ 


= i c — (Xq -T (X\ t. 


/1 2 + (X2 sin/) 2 

By deriving the series in t we obtain the precession rates at J2000.0 

( 72i/f cos/o-|-2(722 /i — ^ 2 ^ 1 ). 

o\=h — H- , - ic 

\ y/(/f + (/2l sin/,,) 2 ) 3 


(/? cot /o + 72, 2 sin 2/o )/? / sin Iq 
CC 1 = i 2 1- , - l c • 


+ sin/o) 2 ) 3 


(42) 


(43) 


( 44 ) 


and 
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For i c = 2.04arc min (Margot et al, 2012) this results in 

5! = —0.00486°/cy and a, = -0.03291°/cy. (45) 

The rotation rate is also slightly modified due to the obliquity. For small i c we get 


I cot I 




(46) 


<p(t) = -M + a- 

/1 2 + (Q sin/) 2 
The resonant rotation rate is consequently 

(3/2) 3 (/i/2i) 2 (3+cos2/ 0 )/2 + (/2 2 /i -M2i)/2i sin2/ 0 + /f/(sin/ 0 ) 2 . 

<Pl = A«0 + ®1 + - , -—- h 

+ {&1 sin/ 0 ) 2 ) 3 

(47) 

and with i c = 2.04arc min this amounts to 6.138506841°/day. The introduced correction is 
not significant when compared to the error of the resonant rotation rate in Eq. 21. 
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